Kosterlitz-Thouless Universality in a Fermionic System 



Shailesh Chandrasekharan and James C. Osborn 
Department of Physics, Box 90305, Duke University, 
Durham, North Carolina 27708, USA 
(Dated: September 21, 2001) 

A new extension of the attractive Hubbard model is constructed to study the critical behavior 
near a finite temperature superconducting phase transition in two dimensions using the recently 
developed meron-cluster algorithm. Unlike previous calculations in the attractive Hubbard model 
which were limited to small lattices, the new algorithm is used to study the critical behavior on 
lattices as large as 128 x 128. These precise results for the first time show that a fermionic system 
can undergo a finite temperature phase transition whose critical behavior is well described by the 
predictions of Kosterlitz and Thouless almost three decades ago. In particular it is confirmed that 
the spatial winding number susceptibility obeys the well known predictions of finite size scaling for 
T < T c and up to logarithmic corrections the pair susceptibility scales as L 2 ~ v at large volumes 
with < n < 0.25 for < T < T c . 
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INTRODUCTION 

A variety of critical phenomena in nature arise due to 
the development of long range correlations or long range 
order at a finite temperature. In the critical region quan- 
tities are characterized by scaling relations that are uni- 
versal and depend only on the symmetries involved. The 
basis for this universality is the renormalization group 
analysis jj| . A number of universal properties have been 
confirmed with great precision by numerically studying 
similar symmetry breaking patterns in different micro- 
scopic models. However, almost all known examples in- 
volve bosonic variables in the form of either classical ||] 
or quantum spins || at the microscopic level. On the 
other hand, most interesting physical systems in con- 
densed matter and high energy physics fundamentally 
involve fermionic particles interacting with each other 
through gauge fields. Unfortunately, calculations con- 
firming the predictions of universality in such systems 
are rather limited and crude. In particular precise calcu- 
lations substantiating the universal critical behavior in 
superconducting materials or strongly interacting mat- 
ter cannot be found. The essential difficulty is that the 
only known approach to study this subject is through 
numerical simulations. In the case of Fermi systems one 
has to be clever and overcome the so called sign problem 
which makes designing algorithms difficult. Even in cases 
where the sign problem can be solved by integrating out 
fermions, the conventional algorithms often suffer from 
critical slowing down. These difficulties have restricted 
the calculations to small system sizes which are not suf- 
ficient to accurately determine the scaling behavior near 
a phase transition. 

The lack of efficient numerical methods available for 
Fermi systems has created some interesting controver- 
sies. It is well known that long range fermionic corre- 
lations can arise at zero temperature and can lead to 



novel universality classes ||. On the other hand it is 
believed that the long range degrees of freedom near a 
finite temperature phase transition should be well de- 
scribed by bosonic degrees of freedom, the reason being 
that fermions acquire a mass proportional to the temper- 
ature T and hence decouple from the critical behavior. 
However, a few years ago this conventional wisdom was 
questioned based on a large N calculation and further 
substantiated by numerical simulations || at smaller N. 
It was suggested that the composite nature of the or- 
der parameter in fermionic systems may cause deviations 
from conventional universality More recently a detailed 
study of the transition revealed that the deviations could 
be related to the artifacts of the large N limit Q. Thus, 
although the controversy for the moment appears settled 
the lesson one learns is that universality classes arising in 
fermionic theories can only be confirmed through precise 
non-perturbative calculations. 

Recently a new method called the meron-cluster algo- 
rithm has emerged as a very efficient alternative to solve 
certain classes of fermionic models [Q, [s|. It is based 
on the well known loop cluster algorithm for quantum 
spin models S and does not suffer from critical slow- 
ing down. Using the concept of a meron-cluster it solves 
the sign problem |hJ. This novel method has lead to 
the first successful determination of the critical behavior 
near a three dimensional Ising transition in a fermionic 
model, again confirming the predictions of universality 
p] , Q. The calculations of M mentioned earlier, had 
reproduced the two dimensional Ising universality class 
from large (even) N fermionic systems. This has now 
been extended to N = 1 fermions using a meron-cluster 
method p3[ . As far as we know, until now no one has 
been able to check the universality arguments with pre- 
cision beyond the Ising universality class starting from 
a microscopic fermionic model. In this letter we present 
the first results from a meron-cluster algorithm which 
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confirms that the critical behavior near a finite temper- 
ature phase transition in a model consisting of fermions 
can be described by the Kosterlitz-Thouless (KT) uni- 
versality class [H. 

The importance of the KT universality class is well 
known. A variety of finite temperature phase transitions 
that arise in many two dimensional condensed matter 
systems involving phenomena like magnetism, supercon- 
ductivity and superfluidity are expected to be described 
by this class. For example the superfluid transition in 
He occurs due to a condensation of Helium atoms which 
in reality are tightly bound objects of fermionic con- 
stituents. Similarly the superconducting transition in 
high T c materials is expected to occur due to a conden- 
sation of electron pairs. Both these transitions are ex- 
pected to follow the predictions of Kosterlitz and Thou- 
less and are related to a U{1) particle number symmetry 
Although KT universality has been studied extensively 
using spin systems, it has been difficult to perform pre- 
cise calculations confirming that a KT universality can 
arise starting from a microscopic fermionic Hamiltonian. 
A simple model that has been studied in this context by 
conventional Monte Carlo methods is the attractive Hub- 
bard model |l5|, Unfortunately the inefficiencies of 
the algorithm have limited the size of the systems that 
could be studied. This coupled with the fact that the 
model has an unusually low critical temperature and ex- 
ponentially diverging correlation length have not yielded 
any quantitative information about the universality class 
of the transition 17 . Given the proliferation of transi- 



tions in nature that are conventionally expected to be 
described by the KT universality class, the fact that at 
a microscopic level most of the systems are always made 
up of fermionic degrees of freedom and the controversy 
discussed above regarding deviations from universality in 
Fermi systems, it is important to find precise calculations 
that either confirm our expectations or contradict them. 



THE MODEL 

The fermionic system we study is a new extension of 
the attractive Hubbard Hamilton operator in two dimen- 
sions and is given by 

H = ^ ^ ^ ^ (Cx,s^ c y,s ~\~ £y,s^Cx.s){\ Sn X y + Tl X y) 



<xy> s=T,| 



l v,i 



4 [S X ■ Sy + J X ■ Jy — J x J y ] 



4 E 



(1) 



Here < xy > represents the nearest neighbor sites of a 
square lattice of size L x L and n xy = n x ^ + n x ^ + 



n y .| + n y> ±. Although the Hamilton operator has many 
terms, it is easy to check that it is invariant under SU (2) 
spin transformations and conserves U(l) fermion num- 
ber. Further, there is an on-site attraction between elec- 
trons of opposite spins like in the attractive Hubbard 
model. As we will see, below a critical temperature trans- 
portation of fermion number through the bulk becomes 
easy, leading to superconductivity (or more appropriately 
superfluidity since the symmetry is not gauged in the 
model) . In higher dimensions this is related to the spon- 
taneous breaking of the U(l) fermion number symme- 
try. In two dimensions, since this is forbidden due to the 
Mermin- Wagner theorem, superconductivity occurs due 
to the KT phenomena. 

The essential feature of the model that makes it useful 
is the fact that its partition function can be written in 
terms of the statistical mechanics of closed loops with 
positive weights (l8) . To see this we start with a discrete 
imaginary time approximation of the partition function 



Z = Tr [exp(-Zf/T)] = Tr 



i=l 



(2) 



where e = l/(TM) and H = H x + H 2 + H 3 + H 4 is a 
convenient reorganization of the interactions present in 
the Hamiltonian. It is then possible to rewrite the path 
integral in terms of closed loops with 



Z = W[b] Sign[6] 

lb] 



(3) 



where [b] is a configuration of bonds linking neighboring 
sites which form closed loop clusters of sites. The magni- 
tude of the Boltzmann weight is W[b] > 0. The fermion 
permutation sign is encoded in the topology of the loops. 
When a configuration contains loops of certain topology, 
referred to as merons, then Sign[6] = 0. If there are no 
meron loops in a configuration then Sign[6] = 2 Nc where 
Nc is the number of loop clusters in the configuration. 
This novel representation of the partition function makes 
the model computationally tractable. More details on the 
meron-cluster formulation can be found in [ fl8| . 

The approximate partition function given in eq. (|^) 
becomes exact in the limit M — > oo. Although this con- 
tinuous time limit can be taken Jl9| , for simplicity the 
results presented here are obtained by fixing M = 20. In 
the study of finite temperature critical behavior a suffi- 
ciently large but finite M is acceptable. This is due to the 
fact that critical dynamics arise from large spatial sizes 
and the discrete nature of the temporal direction only ef- 
fects non-universal quantities. The two dimensional uni- 
versal critical behavior near T r remains unaffected. 
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RESULTS 

The simplest observable relevant to superconductivity 
is the pair susceptibility which we define as 



X 



2T 
~ZV 



l/T 



dt Tr 



-(i/r-t)ff + e -tH - 



(4) 



with p + — ^2 x c x ^c x ^ the pair creation and p~ = (p + Y 
the pair annihilation operators. The susceptibility con- 
tains information about the condensation of electron 
pairs which are the on-site component of Cooper pairs 
in the BCS approach. A formula for the susceptibility in 
terms of loop clusters is easy to construct using the re- 
sults of M . In the present discussed in Q , the 
susceptibility is proportional to the sum over the square 
of the size of certain clusters depending on the number of 
meron clusters in the configuration. In figure [l] we show 
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FIG. 1: Pair susceptibility as a function of L. 

results for the pair susceptibility as a function of spatial 
size L. A state with quasi-long range correlations can 
be seen as a divergence in the susceptibility. Typically, 
for temperatures above the superconducting transition 
temperature T c , the pair susceptibility should reach a 
constant for large enough volumes. This can be seen for 
T = 1.538 and T = 1.429, although in the latter case 
lattices of size L = 128 are necessary to see the satura- 
tion suggesting that the correlation lengths may be on 
the order of 100 lattice units. This dramatic change in 
the correlation length between T — 1.538 and T = 1.429 
is consistent with a Kosterlitz-Thouless prediction that 
it should diverge as exp(Constant/ y/T — T c ) close to the 
critical temperature Q. Interestingly, a fit to a power 
law over the range L < 128 is also poor at T = 1.429. 
Below T c the susceptibility should diverge as 



with the critical exponent r\ starting at 1 /4 at T c and go- 
ing down to as T approaches zero. The exact formula 
has a logarithmic correction to the power law which is 
small and cannot be accurately determined with reason- 
able size lattices so we ignore it here. This continuous 
change in the power is also clearly visible. We find that 
(2 - 77) is 1.767(2) at T = 1.295 and 1.918(3) at T = 0.5. 
The fits to a power law are extremely good for both these 
temperatures over the entire range of L. 
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FIG. 2: Comparison between r/(T) and l/2iv(W 2 } as a func- 
tion of T. 

Although the pair susceptibility alone is sufficient to 
test for the Kosterlitz-Thouless predictions, there could 
always be a lingering doubt whether the power law fits 
would fail when one goes to much larger lattices. For 
example if we only had data for L < 32, a power law 
fit could even work well for T = 1.429 when the error 
bars are sufficiently large. Larger lattices were crucial to 
determine that the power law was not a good fit in that 
case. In order to alleviate such worries we looked at the 
winding number susceptibility which we define as 



(^ 2 ) = <(wy2) 2 + GV2) 2 )/2 



(6) 



where W x (W y ) is the total number of fermions winding 
around the boundary in the x (y) direction. This quan- 
tity is very useful in measuring T c and has been used in 
bosonic systems pi] ]. Although it is a difficult quantity 
to measure with conventional algorithms, it is relatively 
easy in the meron-cluster approach showing the power of 
the new method. Further, we know its finite size scaling 
form to be 
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7T (W 2 ) = 2+ y/A(f) COth( y/A(T) log(L/L (T)) ), 

(7) 

:0@. It can be shown that for T < T c 
00 limit 2irr)(T)(W 2 ) = 1 §1 |1. The 
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fact that (W 2 ) jumps to the universal number 2/n from 
is another well known feature of the KT phenomena 
and can be used to determine T c . Since i](T) can be 
determined from the scaling of the pair susceptibility one 
can combine it with the measurement of (W 2 ) below T c to 
check for consistency in the KT universality class. Figure 
^ compares 1/{2tt{W 2 )) obtained by extrapolating the 
values of the winding number susceptibility to the infinite 
volume limit using the formula (Q) and rj(T) obtained 
from the finite size scaling of the pairing susceptibility, 
as a function of temperature. Clearly the results are in 
excellent agreement with the predictions for a KT phase 
transition. 

DIRECTIONS FOR THE FUTURE 

This work can be extended in several directions. One 
interesting problem in condensed matter physics is to 
understand how disorder effects superconductivity. It is 
predicted that at zero temperature as the disorder is in- 
creased the system undergoes a quantum phase transition 
to an insulating phase p4|. If this is indeed the case, it 
would be interesting to understand the underlying critical 
behavior starting from a fermionic theory. Previous stud- 
ies have used the attractive Hubbard model as the start- 
ing point (25). We suggest that the model studied here 
is perhaps a better alternative since the meron-cluster 
algorithms could turn out to be more efficient in such 
a study. Many other fermionic models with continuous 
symmetries and applications in condensed matter, nu- 
clear and high energy physics can be studied with meron- 
cluster algorithms. A systematic approach has been out- 
lined in fis| . In particular it has also been shown that 
in certain regions of the parameter space some repulsive 
Hubbard-type models with a non-zero chemical potential 
can be solved with meron-cluster techniques. Whether 
there is d-wave superconductivity in these models is an 
open question. A common feature of these new models is 
that they appear to be more complicated. However since 
we have algorithms for them which give a computational 
advantage, they may still produce the clearest results. 
It is a common practice to work with models that are 
the easiest to solve (numerically in this case) to help un- 
derstand complicated physical systems. Once the basic 
phenomena is understood one can then attempt more dif- 
ficult models to explore the robustness and variations in 
the details of the underlying physics. 
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